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Abstract. Based on the behavior of the elastic scattering data, we introduce an almost model-independent 
parametrization for the imaginary part of the scattering amplitude, with the energy and momentum transfer 
dependences inferred on empirical basis and selected by rigorous theorems and bounds from axiomatic 
quantum field theory. The corresponding real part is analytically evaluated by means of dispersion relations, 
allowing connections between particle-particle and particle-antiparticle scattering. Simultaneous fits to 
proton-proton and antiproton-proton experimental data in the forward direction and also including data 
beyond the forward direction, lead to a predictive formalism in both energy and momentum transfer. We 
compare our extrapolations with predictions from some popular models and discuss the applicability of 
the results in the normalization of elastic rates that can be extracted from present and future accelerator 
experiments (Tevatron, RHIC and LHC). 

PACS. 13.85.Dz Elastic scattering - 13.85.-t Hadron-induced high-energy interactions 



X 



1 Introduction 

Elastic hadron-hadron scattering, the simplest hadronic 
collision process, still remains one of the topical theoret- 
ical problems in particle physics at high energies. In the 
absence of a pure QCD description of these large-distances 
scattering states (soft diffraction), empirical analysis based 
on model-independent fits to the physical quantities in- 
volved, play an important role in the extraction of novel 
information, that can contribute with the development of 
useful calculational schemes in the underlying field theory. 

In this context, empirical parametrizations of the scat- 
tering amplitude and fits to the differential cross section 
data have been widely used as a source of model-indepen- 
dent determination of several quantities of interest (the 
inverse problem), such as the profile, the eikonal, the in- 
elastic overlap functions, and, with some additional hy- 
pothesis, even information on form factors (momentum 
transfer space) . These aspects were recently reviewed and 
discussed in pQ , where a list of references to some essen- 
tial results can also be found. However, one aspect of this 
kind of analysis concerns its local description of the ex- 
perimental data, that is, the free parameters are inferred 
from fits to each energy and to each interaction process, 
and therefore the approach has no predictive character. 

In this work we present a novel parametrization for the 
imaginary part of the elastic scattering amplitude with 
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energy and momentum dependences extracted from the 
empirical behavior of the experimental data and selected 
according to some high-energy theorems and bounds from 
axiomatic quantum field theory. The real part of the am- 
plitude is analytically evaluated by means of dispersion re- 
lations, connecting, therefore, particle-particle and particle- 
antiparticle scattering. In this context, the scattering am- 
plitude is expressed as entire functions of the momentum 
transfer and of the logarithm of the energy. Global fits 
to proton-proton (pp) and antiproton-proton (pp) exper- 
imental data in the forward direction (total cross section 
and the p parameter) and, in a second step, also including 
the differential cross sections, lead to a predictive formal- 
ism in the energy and momentum transfer, which is also 
essentially model independent. We present extrapolations 
for values of the energy and momentum transfer above 
those reached in experiments and compare with predic- 
tions from some phenomenological models. We also discuss 
the applicability of the results in the normalization of the 
elastic rates that can be measured in present and future 
accelerator experiments (Fermilab Tevatron, Brookhaven 
RHIC and CERN LHC). 

As will be stressed, this analysis must be seen as a 
first step or attempt toward a formally rigorous model- 
independent description of high-energy elastic hadron scat- 
tering, embodying a predictive character. In this sense we 
shall attempt to discuss and explain, in certain detail, the 
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advantages and disadvantages of the present analysis and 
results. 

The paper is organized as follows. In Sect. |3 we dis- 
cuss the empirical and formal bases of the parametrization 
for the imaginary part of the scattering amplitude and the 
analytical determination of the corresponding real part by 
means of dispersion relations. In Sect.|21we present the fit 
procedures and results, treating firstly only the forward 
scattering and in a second step, including the differential 
cross section data. In Sect. 01 we discuss the physical im- 
plications and applicability of the approach, in the exper- 
imental and phenomenological contexts. The conclusions 
and some final remarks are the contents of Sect. El 



2 Analytical parametrization for the 
scattering amplitude 

The physical quantities that characterize the elastic hadron 
scattering are given in terms of the scattering amplitude 
F, which is expressed as function of two Mandelstam vari- 
ables in the center-of-mass system, usually the energy 
squared s and the momentum transfer squared t = —q 2 - 
We shall base our discussion on the following physical 
quantities the differential cross section, 

^ = T ^\ReF(s ) q 2 ) + ilmF(s,<?)\ 2 , (1) 
dq z lD7rs z 

the total cross section (optical theorem), 
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Fig. 1. Dependence of p from pp and pp elastic scattering (data 
from |5E|1. 

2.1 Parametrization for the imaginary part of the 
amplitude 

2.1.1 Empirical bases 

Let us first investigate some empirical information on the 
differential cross section, in the region of small momentum 
transfer, q 2 < 0.2 GeV 2 . In particular, it is known that at 
q 2 = 0, the data indicate that \p{s)\ < 0.15, as illustrated 
in Fig. 1. Therefore, it is expected that, at small values 
of the momentum transfer, the amplitude F(s,q 2 ) to be 
dominantly imaginary, so that the differential cross section 
in this region can be expressed as 



crtot(s) = 



Im F(s,q 2 = 0) 



(2) 



the p parameter (related with the phase of the amplitude 
in the forward direction), 



P(s) 



ReF{s,q 2 = 0) 
Im F(s,q 2 = 0)' 



(3) 



and the slope of the differential cross section in the forward 
direction, 



(4) 



q 2 =0 



Equations (1) and (2) represent normalizations valid 
in the high-energy region, for example ^/s > 20 GeV [2]. 
We shall return to this point in what follows. 

In this Section we first discuss in certain detail the 
empirical and formal bases that lead to an almost model- 
independent analytical parametrization for the imaginary 
part of the amplitude in terms of both energy and momen- 
tum transfer variables. We then treat the analytical eval- 
uation of the corresponding real part by means of deriva- 
tive dispersion relations and the analytical connections be- 
tween pp and pp scattering. 



da 
dq^ 



1 

16^ 



Im F(s,q 2 ) 



Moreover, in this region, the differential cross section 
data are approximately linear in the logarithm scale, as 
exemplified in Fig. 2, which means that we can express 
the imaginary part of the amplitude as 



Im F{s,q 2 ) 



-N 2 



where a and (3 are real parameters that can depend on 
the energy and reaction considered. 

The point now is to look for possible empirical de- 
pendences for these parameters in terms of the energy s, 
namely analytical expressions for a(s) and f3(s) and that 
is one of the novel aspects of this work. From the above 
two equations and from Eqs. (2) and (4) we have that 
a(s) oc o~tot{s) and (3(s) oc B(s). On the other hand, the 
empirical behavior of <7 4ot (s) and B(s) (near the forward 
direction), displayed in Figs. 3 and 4, respectively, indi- 
cates that in the region of high energies , the empirical 
trends of the data (above y/s « 20 GeV) follow polynomial 
dependences in Ins, of second degree (total cross section) 
and first degree (slope). For these reasons it is reasonable 
to introduce the following empirical parametrizations for 
a(s) and /3(s): 
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Fig. 2. Diffraction peak from pp elastic scattering at y/s = 
52.8 GeV 0. 
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Fig. 3. Experimental information on the pp and pp total cross 
sections from accelerator O and cosmic-ray experiments (see 
|S] for a complete list of references and discussions). 



and 



a(s) = A + B\ns + C In 2 s, 
/9(a) = D + Sins, 



where A, _B, C, £) and E are real constants. We note that 
a dimensionally necessary factor sq, in lns/soi is auto- 
matically absorbed by the other constants. We also note 
that the choice for a(s) is in agreement with the universal 
asymptotic behavior of the total cross sections from the 
analysis developed by the COMPETE Collaboration [§]. 

Now, in the region of medium and large momentum 
transfer, the differential cross section data is characterized 
by the diffractive pattern, as illustrated in Fig. 5. Since we 
have a logarithmic scale, this behavior can be taken into 
account by the standard sum of exponentials in q 2 . 
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Fig. 4. The slope parameter as function of the energy and 
determined in the interval 0.01 < q 2 < 0.20 GeV 2 |7I5) . 
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Fig. 5. Proton-proton differential cross section data at ^/s = 
52.8 GeV 0. 



From the above discussion and aimed to treat both 
pp and pp elastic scattering, we introduce the following 
empirical parametrizations for pp scattering, 



(5) 



with 



a t (s) =A l + B l ln(s) + C, ln 2 (s), 
/3 i (s)=D i + E i \n(s), 



(6) 



and for pp scattering, 



cvi(s)e 



(7) 
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with 

Si (a) ^ A l + B l ln(s) + C t \n 2 (s), (8) 
ft (a) = A + ^M«). 

where i — l,2,...n. In what follows we shall check these 
parametrizations in a formal context. 



2.1.2 Constraints from Axiomatic Quantum Field Theory 



and therefore, in order to not violate the above formal 
results, we must impose the constraint 



i=l 



(9) 



With this condition we also ensure another important 
formal result, namely that if the Froissart-Martin bound 
is reached the p parameter must go to zero logarithmically 

EH 



Even after QCD, unitarity, analyticity, crossing and their 
connections with axiomatic quantum field theory (AQFT) 
still remain a fundamental theoretical framework in the 
investigation of high-energy soft scattering. In this con- 
text, important high-energy theorems and bounds have 
been demonstrated |U)U1 1II12II1.'S| . providing rigorous for- 
mal constraints in the region of asymptotic energies, which 
can not be disregarded in any reliable formalism, mainly 



p(s 



oo oc 



1 

In s 



With parametrizations (5-8) and the constraint (9) we 
have lOn — 1 free parameters, where n is the number of 
exponentials. The novelty in these parametrizations is the 
fact that the energy dependences are already enclosed and 
were inferred from the empirical behavior of the experi- 



related with model-independent approaches. Since parametriz^ental data, being also in agreement with the above high- 



tions (5 - 8) were based exclusively on the behavior of the 
experimental data at fixed (finite) energies, it is necessary 
to check the most important formal asymptotic results. 

Firstly we note that from the optical theorem (2), the 
parametrizations for oti(s) and «i(s) do not violate the 
Froissart-Martin bound, a rigorous prediction of quantum 
field theory which states that 



a tot < cln s, 

where c is a constant. 

Another important result concerns the behavior of the 
difference between particle-particle and antiparticle-parti- 
cle cross sections at the asymptotic regime. In this con- 
text it has been demonstrated by Eden ^01 an d by Grun- 
berg and Truong EH] that if the Froissart-Martin bound is 
reached, the difference between the pp and pp total cross 
sections goes as 



T 1'V 



T 1'P 



Aa 



t' J '' 



< c- 



ln s 



which means that the difference can increases at most as 
In s and even in this case, 



energy theorems. Moreover, the imaginary parts of the 
amplitudes are entire functions of the logarithm of the en- 
ergy, which is an important property in the evaluation of 
the real part, as discussed in what follows. 



2.2 Analytical evaluation of the real part of the 
amplitude 

Connections between real and imaginary parts of the for- 
ward scattering amplitude have been widely investigated 
by means of dispersion relations in both integral and deriva- 
tive forms. In this work we make use of the derivative re- 
lations |17U18| . which are valid in the forward direction 
and for amplitudes belonging to a sub-class of entire func- 
tions of the logarithm of the energy, as it is our case. For 
a recent review and critical analysis on the replacement of 
integral relations by derivative ones see Ref. EH]> where a 
list of references to some outstanding works can also be 
found. 

In the forward direction, the derivative dispersion rela- 
tions for even (+) and odd (— ) amplitudes are expressed 
in terms of a tangent operator and in the case of one sub- 
traction (equal to two subtractions in the even case) they 
are given by [T7IIW??] 



pp 

"tot 1 
pp ) 



as s — > oo (the generalized or revised form of the Pomer- 
anchuk theorem). 

Now, from the optical theorem, Eq. (2) and parametriza- 
tions (5-8) we have 



Aa = J2{(Ai - A,) + (Bi - Bi) In a + (d - &) In 2 a], 



ReF+(s) K 

±±J. = 1- tan 

s s 



Re F_ (a) 



tan 



2 dins 



d 

d In s 



ImF+(s) 



ImF_(s) 



(10) 



(11) 



where K is the subtraction constant. It has also been 
demonstrated by Fischer and Kolaf [20] that at high en- 
ergies the above tangent operator can be replaced by its 
first order expansion, which is the case we are interested 
in. 
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Since besides the forward data we are also aimed to in- 
vestigate differential cross sections, it is necessary to con- 
sider the applicability of the dispersion techniques beyond 
the forward direction. Although several authors make use 
of dispersion relations even for large values of the momen- 
tum transfer, it is important to recall that what is formally 
expected is the validity of the dispersion relations, with a 
finite number of subtractions, inside a region q 2 < q^ax- 
However, the exact expression and/or numerical value of 
Imax depends on the theoretical framework and scattering 
process considered. We shall discuss this subject in some 
detail in Sect. 3.2.2, when applying the formalism to the 
experimental data. Here we only consider a reference to 
this limited interval. 

Based on the above arguments we shall make use of the 
first order derivative dispersion relations, also extended 
beyond the forward region, namely < q 2 < q^ naxl in the 
form 



Rc F + {s,q 2 ) 



K 

s 



A lmF + {s,q 2 ) 



2 d In s 



ReF_(s,g 2 ) tt 



din; 



ImF_(s,9 2 ) 



(12) 



(13) 



Finally, the connections between the hadronic and the 
even/odd amplitudes is established through the usual def- 
initions: 



F pp (s,q 2 ) = F + ( S ,q 2 ) + F_( Sl q 2 ) 
F Pp (s,q 2 ) = F + (s,q 2 ) - F_(s,q 2 ) 



(14) 



This approach is characterized by analytical results 
for both real and imaginary parts of the pp and pp am- 
plitudes. Schematically, from parametrizations (5 - 8) for 
Im F pp / pp (s, q 2 )/ s we obtain Im F + /_ (s, q 2 )/s by inverting 
Eqs. (14). Then the derivative relations (12 - 13) allow to 
evaluate ReF + /_(s,q 2 )/s and by Eqs. (14) we obtain the 
hadronic real parts 



ReF pp (s,q 2 ) K 



+ E{f His) - cciismsW] e'^W 



~ [«i(s)< 



-Ms)q 2 



«i(s)e" 



-ft(s)? 2 



}. 



ReF Pp {s,q 2 ) _ K 



n 



-ft(s) 9 2 



where the primes denote differentiation with respect to 
Ins (Eqs. (6) and (8)). With this we have analytical ex- 
pressions for the pp and pp differential cross sections: 



da 



pp/pp _ 1 
~ 16tt 



ReFpp/ppis^q 2 ) Im F pp/pp (s,q 2 ) 



dq 



It should be noted that exact analyticity and crossing 
properties demand symmetric variables, namely the labo- 
ratory energy E for q 2 — and the variable (s — u)/4m, 
where u is the Mandelstam variable, for q 2 > 2 . How- 
ever, since E depends linearly on s and, as will be dis- 
cussed in Sect. 3.2, we shall consider the applicability of 
the formalism mainly in limited regions of the momentum 
transfer and only above y/s — 20 GeV, the use of s as 
variable does not introduce essential changes in the above 
formulas. 

Taking into account the subtraction constant K, the 
constraint (9) and parametrizations (5-8) we eventually 
have lOn fit parameters in the case of n exponential terms. 
This completes the analytical construction of the formal- 
ism, characterized by its empirical basis, essentially model- 
independent parametrizations, agreement with high-energy 
theorems and amplitudes belonging to the class of entire 
functions in the logarithm of the energy. In the next sec- 
tion we determine the free parameters involved through 
fits to pp and pp elastic scattering data. 



3 Experimental data, fitting and results 
3.1 Experimental Data 

The most important empirical input in our parametriza- 
tion is the energy dependence enclosed in the expressions 
of a(s) and /3(s), Eqs. (6) and (8), respectively. Since it 
characterizes the region where the total cross section in- 
creases with the energy, we shall consider here only the 
experimental data available above s/s = 20 GeV from pp 
and pp scattering. We note that this necessary threshold 
puts limitations in extensions of the formalism to other re- 
actions, such as ir^p, K^p, etc., due to the small number 
of experimental data available. 

For the forward data on atot and p, we use the Particle 
Data Group archives [3], to which we added the value of 
p and atot at 1.8 TeV obtained by the E811 Collaboration 
0. The statistical and systematic errors were added in 
quadrature. We did not include the cosmic-ray informa- 
tion on pp total cross sections due to the model depen- 
dences involved [HI- 

The differential cross section data include the optical 
point, 



da(s, q 2 
dq 2 



0-*ot(l 



o 2 ) 



16tt 



(15) 



and the data above the Coulomb-nuclear interference re- 
gion, namely q 2 > 0.01 GeV 2 . The data include 12 sets 
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form pp scattering, at y/s = 23.5, 27.4, 30.7, 44.7, 52.8, 
and 62.5 GeV and from pp scattering, at y/s — 31, 53, 61, 
546, 630 and 1800 GeV. The pp data at 27.4 GeV, cover- 
ing the region 5.5 < q 2 < 14 GeV 2 , are from [23- The pp 
differential cross section data at 1.8 TeV include those ob- 
tained by the E710 Collaboration |S] (0.045 < q 2 < 0.627 
GeV 2 ) and by the CDF Collaboration [221 (°- 035 < I 2 < 
0.285 GeV 2 ). In this case we used two optical points with 
the values of a tot and p from references ^ (E811 Collab- 
oration) and [Hj (E710 Collaboration). The complete list 
of references to the other data sets can be found in 1 
(references [26, 28-31]). In all these sets the experimental 
errors correspond to the statistical ones. 

We note that we have used all the experimental data 
referred to before, that is, we did not performed any kind 
of data selection in the above standard ensemble. 



3.2 Fitting and results 

As recalled in Sec. 2.2, the applicability of the disper- 
sion relations, outside the forward direction, depends on 
the maximum value of the momentum transfer considered. 
For this reason, we shall treat separately fits to only for- 
ward quantities, atot and p, and simultaneous fits to these 
quantities plus the differential cross section data. We first 
present the fits to the forward data and next discuss the 
applicability of the dispersion techniques beyond the for- 
ward direction. 



3.2.1 Fits to the forward scattering data 

Making use of the formalism described in Sect. 2, we per- 
formed simultaneous fits to <j t ot and p data, above 20 GeV, 
from pp and pp scattering. Since we are treating here only 
forward data (q 2 = 0), the sole parameters involved are 
those associated with aj(s) and ati(s) in Eqs. (6) and (8). 

The fits were performed through the CERN-Minuit 
code, with the estimated errors in the free parameters cor- 
responding to an increase of the % 2 by one unit. For this 
ensemble of data good statistical results were obtained 
with only one exponential factor (n = 1 in Eqs. (5) and 
(7)) and the best fit indicated % 2 /DOF =1.07 for 83 de- 
grees of freedom. The constraint (9) in this case reduces 
to Ci = Ci and the fit indicated a value of the subtrac- 
tion constant compatible with zero. The numerical results 
are displayed in Table 1 and the corresponding curves, to- 
gether with the experimental data analyzed, are shown in 
Figure 6. 

These results will be discussed in detail in Sect. 4, but 
we note here the good quality of the fit in terms of the 
\ 2 /DOF and also the small number of free parameters 
involved: 5. We also note a crossing in the total cross sec- 
tions, with crfo t becoming higher than af^ t above y/s ~ 
100 GeV, and a similar effect is predicted for p(s). As we 
recalled in Sect. 2.1, these behaviors do not violate any 
high-energy theorem on elastic hadron scattering. How- 
ever, the result for pP p (s) is below the experimental data 
available at the highest energies. We shall discuss this ef- 
fect in Sec. 4.2. 
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Fig. 6. Fit results to the forward data, atot and p, from pp 
and pp scattering. 



3.2.2 Fits beyond the forward direction 

We now consider simultaneous fits to a to t{s), p(s) and 
d<T(s,q 2 )/dq 2 , from pp and pp scattering. As recalled be- 
fore, although dispersion relations have been used even 
in the region of medium and large momentum transfer 
(see, for example, [18 , 23 , 24 25 . 25]), an important point 
concerns the exact region in the q 2 variable inside which 
dispersion relations hold. In what follows we first review 
some formal results involved which show us that, in the 
case of pp and pp scattering (and nucleon-nucleon in gen- 
eral), the situation is not simple or neat. Based on these 
results we shall infer a reasonable strategy, not proved to 
be wrong, that will allow us to develop simultaneous fit 
procedures including the differential cross section data. 
Note that, since we are treating with a sub-class of en- 
tire functions in the logarithm of the energy, the discus- 
sion that follows applies equally well to both integral and 
derivative dispersion relations |19) . 



• Analyticity in q 2 

Dispersion relations are connected with the unitarity 
and analyticity properties of the amplitude. Recently, the 
axiomatic approach to high-energy hadron scattering, as 
well as the rigorous analyticity-unitarity program have 
been nicely reviewed in the excellent papers by Vernov 
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Table 1. Results of the simultaneous fits to a to t and p from pp and pp scattering. All the parameters are in GeV 2 , x 2 /DOF 
= 1.07 for 83 degrees of freedom and C\ =G\. 





pp scattering 




pp scattering 


Ai 


121.9 ± 2.7 


Ai 


140.8 ± 3.6 


Si 


-9.82 ± 0.72 


Si 


-11.78 ± 0.86 


Ci 


1.036 ± 0.049 


C x 


1.036 ± 0.049 



and Mnatsakanova an d Martin ^3] > where a complete 
list of references and credits to outstanding results and au- 
thors can be found. For this reason, based on these works, 
we shall only summarize and quote here the results, we 
understand, give an updated view on the (^-interval in- 
side which dispersion relations hold. The results are the 
followings. 

1. For meson-meson and meson- nucleon scattering, rig- 
orous formalism based on local field theory, allows to 
prove that dispersion relations are valid in finite inter- 
vals of the momentum transfer. For example, g^az = 
32m 2 /3 w 0.2 GeV 2 in the case of ttN scattering [12] 
and q 2 lax = 28m 2 m 0.55 GeV 2 in the case of hit scat- 
tering [T5| . 

2. Similar intervals can be deduced for other processes, 
like 7 + N -> 7 + ATW, 7 + N -> tt + ATM, e + N -> 
e + 7r + iVW [T3]. 

3. In the formal analyticity-unitarity context there seems 
to be no results for nucleon-nucleon scattering. How- 
ever, a limit q^ax = m 1/^ ~ 0.005 GeV 2 can be in- 
ferred from perturbation theory |13| . 

4. The reason why elastic pp and pp amplitudes "lack 
the usual analytical properties is that the cut in the 
complex s plane starts from sq — 4m 2 (due to virtual 
annihilation process), while the physical region of pp 
scattering starts from si = 4M 2 " Here m corre- 
sponds to the pion mass and M to the proton mass. 

5. In the context of the double-dispersion representation 
by Mandelstam [2Z|, the domain in the q 2 variable, 
inside which dispersion relations hold for a process 
m + m — > m + m, extends up to q 2 nax = 9m 2 |28| . 
Although the original approach treated only pion-pion 
scattering, one point to stress is the fact that for all 
mass cases this representation was never proved nor 
disproved in the contexts of the axiomatic field theory 
or perturbation theory |13j . 

We also recall that fixed-g 2 dispersion relations for 
nucleon-nucleon scattering have been used by Kroll and 
co-authors [2*4ll2*5] and in particular, in Ref. j2U, pp and 
pp scattering were investigated through dispersion rela- 
tions in the region < q 2 < 3 GeV 2 . However, there is no 
reference to a formal numerical value for <£n OX . 

These are the results we have found and compiled on 
the applicability of dispersion relations beyond the for- 
ward direction. We understand that the quoted bound 
from perturbation theory seems unreliable to be consid- 
ered in the case of a soft processes like elastic pp and 
pp scattering. A second aspect concerns the fact that, in 
the context of the axiomatic field theory, the Mandesltam 



representation for all mass cases was never proved to be 
incorrect (or correct either). If the representation can be 
extended to the pp and pp case, the analyticity domain 
could cover the region up to q^y iax — 9m 2 w 8 GeV 2 . 

• Strategies and fits 

Based on the above information, we understand that 
it can be instructive to perform tests on distinct values 
for q^ax and investigate the consequences in the descrip- 
tion of the bulk of the experimental data on atot(s), p(s) 
and da(s,q 2 )/dq 2 , from pp and pp scattering. Since the 
typical mass scale in the hadronic scattering is the proton 
mass (which is also expected to represent an interface be- 
tween soft and semihard processes), it may be reasonable 
and perhaps even conservative, to consider some bounds 
q^ax inside the region 1-2 GeV 2 . Moreover, it seems 
also important to address the practical applicability of 
the dispersion approach at medium and large values of 
the momentum transfer by taking into account all the dif- 
ferential cross section data available, namely, 9^ ax = 14 
GeV 2 . Despite the lack of a formal justification for this 
extreme case, we understand that it may also be useful 
to get some additional information on the regions where 
dispersion relations work, even if only in an strictly phe- 
nomenological context. It is important to stress that the 
strategy to consider different values for g^ ax is only an 
ansatz and that the main point in favor of this hypothesis 
is the fact that there is no formal proof against it, or in 
other words, we understand that it should not constitute 
a serious formal drawback. 

Based on the above discussion, we shall consider four 
variants for the fits by selecting differential cross section 
data up to g 2 nax = 1.0, 1.5, 2.0 and 14 GeV 2 . As before, the 
fits were performed through the Minuit code. For these en- 
semble of data, independently of the value considered for 
<7^ ax , the best results demanded three exponential terms 
in the imaginary part of the amplitude and therefore 30 
free parameters to be fitted. The constraint (9) was taken 
into account by defining C\ = C\ + Ci + C3 — Ci — C3. 

The x 2 information on each of the four variants consid- 
ered is displayed in Table 2. We note that the x 2 /DOF lies 
in the interval 2.5 - 3.0 for a number of degrees of freedom 
equal or greater than 923. It is important to mention that 
these values are typical of global fits to the experimen- 
tal data on atot(s), p(s) and do(s, q 2 )/dq 2 , from pp and 
pp scattering |29|. The "large" values are consequences of 
several points in the differential cross sections that lies 
outside a normal distribution, as well as different normal- 
izations from different experiments in distinct kinematic 
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Table 2. Statistical information on the fit results to a to t, p and da/dq 2 data from pp and pp scattering in different intervals of 
the momentum transfer variable. 



?L GeV 2 


NDOF 


X 2 /DOF 


f.O 


923 


2.476 


1.5 


1003 
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2.0 
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14.0 (all data) 
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2.829 
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Fig. 7. Differential cross sections from global fits to pp and pp 
data with q^ax = 2 GeV 2 . Curves and data were multiplied by 
factors of 10 ±2 . 
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Fig. 8. Total cross section and p from global fits to 
data with differential cross section data up to g^ax : 



pp and pp 
= 2 GeV 2 . 



intervals. As commented before we did not perform any 
kind of data selection. 

In fact, despite the large values of the x 2 /DOF, the vi- 
sual description of the experimental data is good in all the 
cases investigated. In particular we display here the results 
for g^ ax = 2 GeV 2 , which we understand can be consid- 
ered a conservative case (in agreement with the expected 
analyticity interval in terms of the momentum transfer) 
and for q^ ax = 14 GeV 2 . The values of the free parameters 
in both cases are shown in Table 3 and the corresponding 
curves together with the experimental data analyzed in 
Figs. 7 and 8 (g^ ax = 2 GeV 2 ) and 9 and 10 (<?^ ax = 14 
GeV 2 ). We see that the description of all the differential 
cross section data is quite good, even in the extreme case 
q 2 max = 14 GeV 2 . 



4 Discussion 

In this Section we first summarize the main results we 
have obtained and then proceed with a discussion on their 
physical implications and their applicabilities in the exper- 
imental and phenomenological contexts. 

By means of a novel essentially model-independent an- 
alytical parametrization for the scattering amplitude and 
fits to physical quantities that characterize the elastic pp 
and pp scattering, we have developed a predictive formal- 
ism in the variables s and q 2 , that has only empirical 
and formal bases. The approach is intended for the high- 
energy region, specifically above -y/s — 20 GeV (in order 
to guarantee the empirical energy dependences). We first 
considered global fits to only the forward data, a to t and 
p, for which dispersion relations can be formally applied. 



Avila, Campos, Menon, Montanha: Phenomenological analysis connecting pp and pp elastic scattering 



9 



Table 3. Results of the simultaneous fits to a to t, p and da/dq 2 from pp and pp scattering, with differential cross section data up 
to q 2 max = 2.0 GeV 2 and q^ ax = 14 GeV 2 , for which K = 49. 7± 1.7 and K = -0.10534 0.0048, respectively. All the parameters 
are in GeV" 2 and d = C\ 4 C 2 4 C 3 - C 2 - C 3 . 



pp scattering pp scattering 





q 2 = 2.0 GeV 2 


qLax = 14 GeV 2 

M f MUX 




q 2 = 2.0 GeV 2 

triad- 


qLax = 14 GeV 2 




91.13 ± 0.28 


109.70 ± 0.28 


Ai 


119.61 ± 0.44 


112.28 ± 0.44 






i r v>q 4- n rnq 

lu.jiy zn u.uoy 


Jz>i 


9 4SR 4 D fl7^! 
i.'iou zn u.u ( o 




Ci 


constrained 


constrained 


d 


-0.0174 ± 0.0038 


-0.1673 ± 0.0039 


Dl 


-7.79 ± 0.33 


-8.91 ± 0.32 




3.134 ± 0.067 


3.170 ± 0.069 


Si 


2.908 ± 0.051 


3.045 ± 0.050 


E x 


0.4884 ± 0.0078 


0.4860 ± 0.0082 


A 2 


16.82 ± 0.22 


-4.06 ± 0.23 


A 2 


14.51 ± 0.15 


10.23 ± 0.15 


#2 


7.071 ± 0.030 


11.387 ± 0.030 


B 2 


-6.730 ± 0.024 


-6.756 ± 0.027 


c 2 


0.3027 ± 0.0047 


0.0952 ± 0.0047 


c 2 


0.9035 ± 0.0027 


0.9613 ± 0.0029 


D 2 


1.647 ± 0.014 


1.290 ± 0.014 


D 2 


-1.549 ± 0.011 


-1.476 ± 0.013 


E 2 


0.4646 ± 0.0022 


0.5097 ± 0.0023 


E 2 


0.5521 ± 0.0011 


0.5645 ± 0.0012 


A 3 


0.2582 ± 0.0049 


1.0554 ± 0.0077 


A z 


-17.160 ± 0.055 


-8.148 ± 0.031 


B 3 


-0.09894 ± 0.00074 


-0.3607 ± 0.0013 


B z 


2.6083 ± 0.0060 


1.2313 ± 0.0040 


C z 


0.5921E-02 ± 0.0070E-02 


0.02372 ± 0.00012 


Cz 


-0.11298 ± 0.00038 


-0.05572 ± 0.00025 


Dz 


0.303 ± 0.019 


0.6454 ± 0.0081 


Dz 


1.5672 ± 0.0084 


0.9272 ± 0.0072 


Ez 





0.0176 ± 0.0011 


E 3 





0.03859 ± 0.0011 




Fig. 9. Differential cross sections from global fits to pp and pp 
data with q 2 ^^ = 14 GeV 2 (all differential cross section data). 
Curves and data were multiplied by factors of 10 ±2 . 
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Fig. 10. Total cross section and p from global fits to pp and 
pp data with q^a* = 14 GeV 2 (all the differential cross section 
data). 
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We then included the differential cross section data and 
discussed strategies for the use of the dispersion relations, 
namely fits in different intervals in the momentum transfer 
variable: q 2 max = 1.0, 1.5, 2.0 and 14 GeV 2 (all data). The 
main point is the fact that there is no formal proof against 
these assumptions. Although we have displayed here only 
the results for q 2 = 0, q 2 max = 2 and q 2 max =14 GeV 2 , in 
all the cases investigated we have obtained good descrip- 
tions of the experimental data analyzed. As commented 
before, we consider the results with g^ aa . =14 GeV 2 as an 
illustrative example on the practical applicability of the 
dispersion relations at medium and large values of the mo- 
mentum transfer. However, a striking feature is the high 
quality of the data description reached in this case, as 
shown in Figs. (9) and (10). 

In what follows, we discuss the applicability of the 
physical results in the experimental and phenomenological 
contexts. In the former case we shall consider processes 
that are being investigated or planned to be treated in 
accelerator experiments, referring to the following three 
cases: (1) pp scattering at y/s = 200 GeV, that was in- 
vestigated and might yet be investigated by the pp2pp 
Collaboration at the Brookhaven RHIC; (2) pp scatter- 
ing at y/s — 1.96 TeV, that are being analyzed by the 
DZcro Collaboration at the Fermilab Tevatron (RUN II); 
(3) pp scattering at y/s = 14 TeV, planned to be investi- 
gated by the TOTEM Collaboration at the CERN LHC. 
In the phenomenological context we shall make reference 
to some popular models which are, at the same time, 
representatives of different pictures to high-energy soft 
diffraction. To this end we will limit the discussion here 
only to the models by Desgrolard, Giffon and Predazzi 
(DGP) H3, Bourrely, Soffer and WU (BSW) gOj, Don- 
nachie and Landshoff (DL) Block, Gregores, Halzen 
and Pancheri (BGHP) [32] and the Odderon concept, in- 
troduced by Lukaszuk and Nicolescu |33| . 

We shall focus our discussion on the predictions ob- 
tained for q 2 — (forward data only), q^ ax — 2 GeV 2 
and g^ OI = 14 GeV 2 . Also, we treat separately the re- 
sults for the total cross section, the p parameter and the 
differential cross section, obtained from the above three 
variants of the fit procedure. 

4.1 Total cross section 

In the three fit variants, the results indicate a crossing with 
a PP t becoming greather than <J PP t . However, the crossing 
point depends on the g 2 -interval of the differential cross 
section data considered in the fit: y/s ps 100 GeV for q 2 = 
0, « 500 GeV for q 2 max = 2 GeV 2 and ps 3 TeV for q 2 max = 
14 GeV 2 as shown in Figs. 6, 8 and 10, respectively. That 
is, asymptotically, the difference Aa = a\ p t — a p t p t does not 
goes to zero and as recalled in Sec. 2.1, this behavior is 
not in disagreement with formal results obtained in the 
context of axiomatic field theory. 

Among the models quoted, this result suggest a dom- 
inant contribution of the Odderon at the highest en- 
ergies. We stress that, the only information we have in- 
troduced in our paramctrization was the empirical fact 



that the pp and pp total cross sections increase at most 
as In 2 s (parametrizations for o^(s) and 5j(s) in Eqs. (6) 
and (8)) and that the difference may increase at most as 
Ins (constraint (9)). 

In particular, at y/s = 1.80 TeV, the experimental re- 
sults for a PP t are characterized by the well known discrep- 
ancies between the values reported by the E811 and E710 
Collaborations 0flS] and that reported by the CDF Col- 
laboration In this respect, except for the forward fit 
result for afg t , which lies between the discrepant points 
(Fig. 6), the predictions including the differential cross 
section data favor the E811/E710 results (Figs. 8 and 10). 

In Table 4 we present our numerical predictions for the 
total cross sections in the case of the experiments referred 
to before. For pp scattering at 14 TeV (LHC) our results 
with = 2 GeV 2 and q^ax — 14 GeV 2 are in agree- 

ment, respectively, with the predictions from the BGHP 
model (cr^ = 108.0 ± 3.4 mb) and from the BSW 
model {a tot = 103.5 mb) [23j. However it should be noted 
that these models do not distinguish pp and pp scattering 
at asymptotic energies. The Table also contains the results 
for p(s), to be discussed in what follows. 

4.2 The p parameter 

As a consequence of the connections between real and 
imaginary parts of the amplitudes via dispersion relations, 
similar effects appear in our results for p(s), as shown in 
Figs. 6, 8 and 10: p pp {s) becomes greater than pP p {s) above 
yfs « 80 GeV (q 2 = 0), « 200 GeV {q 2 max = 2 GeV 2 ) and 
ss 2 TeV (q^ax = 14 GeV 2 ). In all the cases the constraint 
(9) assures the asymptotic behavior as 1/lns for both pp 
and pp scattering. 

As in the case of the total cross section, these re- 
sults are in agreement with the Odderon dominance at 
the highest energies. A crossing in p(s) with p pp (s) be- 
coming greather than pP p (s) is also predicted in one of 
the versions of the DGP model (2H! and in the analysis of 
Rcf . 6 , which includes cosmic-ray information on a PP t and 
a model-dependent parametrization with Odderon contri- 
bution. 

Differently from the results for the total cross sections, 
we note here some distinct characteristics between the pre- 
dictions for p(s) obtained with only the forward data (Fig. 
6) and those including the differential cross section data 
up to q 2 max = 2 GeV 2 (Fig. 8) and q 2 max = 14 GeV 2 (Fig. 
10). In the former case the curve for p pp (s) lies below the 
highest pp data, what does not occur when the differential 
cross section data is included. We have realized that this 
effect (Fig. 6 and partially in Fig. 8) is due to the large 
error bars of the experimental data at 1.8 TeV and also to 
the small number of p data from pp scattering above 20 
GeV. In fact, at y/s — 1.8 TeV, the experimental values 
are: p E811 = 0.132 ± 0.056 @J and p E 7io = 0.140 ± 0.069 

corresponding to relative errors of 42% and 49%, re- 
spectively. For example, if we use the same central values 
and reduce the errors to 10 %, the same fit leads to a 
curve that pass through the central values. However this 
is only a technical information that certainly has nothing 
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Table 4. Predictions for atot(s) and p(s) from fits to only forward data (q 2 = 0) and including the differential cross section 
data (q 2 max = 2 GeV 2 and q 2 max = 14 GeV 2 ). 





Process 







2 

*7max 


2 GeV 2 


2 

*7max 


14 GeV 2 






atot (mb) 


P 


a to t (mb) 


P 


a t ot (mb) 


P 


pp - 


y/s = 200 GeV 


52.27 


0.1532 


51.32 


0.1439 


51.12 


0.1065 


pp - 


Vs = 1.96 TeV 


78.05 


0.1114 


75.74 


0.1124 


75.24 


0.1343 


pp - 


yfs = 14 TeV 


121.6 


0.1964 


107.5 


0.1321 


105.4 


0.1365 



to do with a physical result in the context of our analysis. 
Experimentally it is known that, as the energy increases, 
it is very difficult to reach the Coulomb-nuclear interfer- 
ence region, from which the p parameter is extracted 
Therefore it is not expected an improvement in these ex- 
perimental values, unless some novel technique could be 
developed. In this respect, the above effect, at the highest 
energies, can not be eliminated in the present formulation 
and fit procedure, constituting, therefore, a drawback in 
our analysis, when only forward data is considered. How- 
ever, comparison of Figs. 6, 8 and 10 shows an interesting 
effect: the quality of the visual description of the p data 
at the highest energies is improved with the addition of 
the differential cross section information. 

The numerical predictions for p(s), in the case of the 
experiments referred to before, and from the three fit vari- 
ants (q 2 — 0, q 2 nax — 2 GeV 2 and g 2 aaa . = 14 GeV 2 ) are 
displayed in Table 4. 

4.3 Differential cross section 

As we have shown, the descriptions of the pp and pp differ- 
ential cross section data analyzed are quite good for 
= 2 GeV 2 (Fig. 7) and even in the case of q 2 max = 14 GeV 2 
(Fig. 9). In this subsection we discuss the applicability of 
these results in the experimental and phenomenological 
contexts. 

To the extend that our analysis can be considered 
model-independent and predictive, it may be instructive 
to detail the results for the experiments referred to in the 
beginning of this Section. Our predictions for these pro- 
cesses, from the global fit including the differential cross 
section data up to q^ax = 2 GeV, are shown in Fig. 11 
and the corresponding numerical results, in the region 
- 2 GeV 2 , are displayed in Table 5 for some values of the 
momentum transfer. From Fig. 11 we note the presence of 
a dip at q 2 w 1.2 GeV 2 for pp scattering at 200 GeV and 
that the diffraction pattern becomes a shoulder at higher 
energies for both pp (y/s = 1.96 TeV) and pp scattering 
(•y/s = 14 TeV); we also note the shirinkage of the diffrac- 
tion peak as the energy increases. 

4.3.1 Experimental aspects 

As examples of practical use of these predictions in the 
experimental context, let us discuss the recent determina- 
tions of the slope parameter from elastic rates measured 





pp,Vs=14.0 TeV 




pbarp,Vs=1 .96 TeV 




- - pp,Vs=200.0 GeV 











0.0 0.5 1.0 1.5 2.0 

Fig. 11. Predictions for the differential cross sections at the 
RHIC, Tevatron and LHC energies from fits including the dif- 
ferential cross section data up to q^a^ = 2 GeV 2 (Table 5). 
The upper and lower curves were multiplied by 10 3 and 10 -3 , 
respectively. 



by the pp2pp Collaboration {pp scattering at s/s = 200 
GeV) at the RHIC [35] and the preliminary results ob- 
tained by the DZero Collaboration (pp scattering at \fs = 
1.96 TeV) at the Tevatron Operationally, the differ- 
ential cross section is expressed by 



da _ 1 dN 
dq 2 C dq 2 ' 

where dN/dq 2 is the rate of the clastic interactions and 
C the machine luminosity. Due to uncertainties in the de- 
termination of £ the above quoted experiments have ex- 
tracted only the slope of the elastic rates, that is, the 
corresponding differential cross section could not yet be 
determined. In what follows we present our results for the 
corresponding slopes and discuss ways to contribute with 
a possible reasonable normalization of the elastic rates. 

• pp at «/s = 200 GeV 

In the case of the pp2pp experiment, the slope B was 
obtained from the elastic rates measured in the q 2 range 
0.010 < q 2 < 0.019 GeV 2 . The corresponding amplitude 
has contributions from the Coulomb amplitude, nuclear 
amplitude and the interference between them, and it is 
parametrized by 
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Table 5. Predictions for the differential cross sections in mbGeV~ 2 at the RHIC, Tevatron and LHC energies, from global fits 
including the differential cross section data up to g^ a:c = 2.0 GeV 2 . 



q 2 GeV 2 


pp - y/s = 200 GeV 


pp - y/s= 1.96 TeV 


pp - y/s = 14 TeV 


0.00 


136.24 


298.18 


602.87 


0.01 


117.66 


251.98 


470.45 


0.05 


66.649 


129.03 


190.54 


0.10 


33.576 


56.399 


65.903 


0.15 


17.144 


24.839 


23.222 


0.20 


8.8167 


10.982 


8.2384 


0.25 


4.5585 


4.8519 


2.9339 


0.30 


2.3699 


2.1307 


1.0537 


0.35 


1.2400 


0.92503 


0.39265 


0.40 


0.65358 


0.39539 


0.16557 


0.45 


0.34731 


0.16687 


0.090879 


0.50 


0.18607 


0.071284 


0.067705 


0.55 


0.10042 


0.033293 


0.060450 


0.60 


0.054478 


0.019393 


0.057098 


0.65 


0.029605 


0.014952 


0.054020 


0.70 


0.016034 


0.013790 


0.050435 


0.75 


0.0085958 


0.013451 


0.046423 


0.80 


0.0045185 


0.013078 


0.042242 


0.85 


0.0022989 


0.012456 


0.038117 


0.90 


0.0011103 


0.011606 


0.034201 


0.95 


0.00049312 


0.010614 


0.030578 


1.00 


0.00018982 


0.0095647 


0.027284 


1.05 


5.5411E-05 


0.0085216 


0.024324 


1.10 


8.5222E-06 


0.0075277 


0.021686 


1.15 


3.9878E-06 


0.0066078 


0.019345 


1.20 


1.7294E-05 


0.0057734 


0.017276 


1.25 


3.5677E-05 


0.0050275 


0.015449 


1.30 


5.2974E-05 


0.0043677 


0.013838 


1.35 


6.6641E-05 


0.0037883 


0.012417 


1.40 


7.6041E-05 


0.0032825 


0.011164 


1.45 


8.1462E-05 


0.0028425 


0.010058 


1.50 


8.3576E-05 


0.0024608 


0.0090810 


1.55 


8.3151E-05 


0.0021304 


0.0082164 


1.60 


8.0907E-05 


0.0018448 


0.0074507 


1.65 


7.7458E-05 


0.0015981 


0.0067714 


1.70 


7.3291E-05 


0.0013851 


0.0061679 


1.75 


6.8779E-05 


0.0012013 


0.0056310 


1.80 


6.4189E-05 


0.0010426 


0.0051523 


1.85 


5.9707E-05 


0.00090562 


0.0047250 


1.90 


5.5453E-05 


0.00078734 


0.0043428 


1.95 


5.1497E-05 


0.00068517 


0.0040002 


2.00 


4.7875E-05 


0.00059687 


0.0036928 



B = 16.3 ± 1.6 (stat.) ± 0.9 (syst.) GeV" 2 

Adding in quadratures the error reads ± 1.8 GeV -2 . 

From Fig. 4 we can see that this experimental value of 
the slope is above the general trend of the others measure- 
ments, even from pp scattering. This effect is due to the 
small values of the momentum transfer in which the mea- 
surement has been performed, namely lower than those in 
the other experiments and also because the interval is in 
the limit of the Coulomb-nuclear interference region (q 2 w 
0.01 GeV 2 ). Since the pp data we have analyzed cover the 



— - 4n(fic) 2 ( + 1 + p2 j 2 e - s l*l 
dt -^ hc > { t ) + l 6 7T(?ic) 2 tot 

The fit parameters are the slope B and a normaliza- 
tion constant (elastic rates). The input values for a to t and 
p used by the authors were 51.6 mb (obtained from the 
Donnachic-Landshoff model) and 0.13 (fit by the UA4/2 
Collaboration), respectively. The resulting slope parame- 
ter was 
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region only up to yfs — 62.5 GeV and above q 2 = 0.01 
GeV 2 (except for the optical point), it is an important 
test to check our predictions for the above quantity. 

To this end, from the fit with q 2 lax = 2 GeV 2 and based 
on the experimental procedure |35|. we have generated 19 
differential cross section points, with estimated error of 
1%, in the region 0.010 < q 2 < 0.019 GeV 2 and fitted 
the points with an exponential form in the momentum 
transfer 



da 
dj 2 



(16) 



as shown in Fig. 12. With this procedure we have obtained 



A = 136.0 ± 1.7 mbGeV , 
B = 14.46 ±0.84 GeV -2 , 

with x 2 /DOF = 4.8 x 10~ 5 for 17 degrees of freedom. 
Therefore, our result for the slope is in agreement with 
the experimental value, lying inside the lower error bar in 
the case that statistical and systematic errors are added 
in quadrature. The relative error in respect the central 
value is 11%. Moreover, the input value used by the pp2pp 
Collaboration, a\^ t = 51.6 mb, is also in agreement with 
our predictions for the cross section, as shown in Table 4, 
namely 51.32 mb. We also note that, although our results 
indicate p = 0.1439, which is higher than the input p = 
0.13, this difference has no practical effect in the nuclear 
contribution in Eq. (16), since this parameter appears in 
the form 1 + p 2 . However that is not the case for the total 
cross section which has a quadratic contribution: a 2 ot . We 
understand that these results corroborates the accuracy 
of our predictions and in this sense, the above value we 
obtained for the parameter A could be used as a suitable 
normalization factor in the estimation of the correspond- 
ing differential cross section. 



• ppatyfs = 1.96 TeV 

Now let us discuss the recent measurements (even if 
preliminary) of the elastic rates performed by the DZero 
Collaboration, from pp scattering at yfs = 1.96 TeV |36| . 
In this case, the rate of elastic collisions has been mea- 
sured at medium values of the momentum transfer, in the 
interval 0.96 < q 2 < 1.31 GeV 2 [SZj. As an illustration, 
and for further discussion, our prediction for the differen- 
tial cross section at yfs = 1.96 TeV and yfs — 1.80 TeV, 
with q 2 lax = 2 GeV 2 , are shown in Fig. 13 together with 
the experimental data, obtained by the E710 and CDF 
Collaborations, at yfs =1.80 GeV. 

In principle, the elastic rates at 1.96 TeV could be com- 
pared with the differential cross section data at 1.80 TeV, 
allowing a kind of normalization. However, from Fig. 13, 
we see that the E710 data cover the region only up to q 2 
= 0.627 GeV 2 and the main problem is the fact that in 
the gap between this last point and the first DZero point 
(q 2 « 0.96 GeV 2 ) it is expected the presence of a dip or a 
shoulder, implying, in any case, in a change of curvature. 
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What is worst, the E710 point at q 2 = 0.627 GeV 2 has 
a large error bar (not shown in the figure, but taken into 
account in all the fits), turning out very difficult, in our 
oppinion, any attempt to perform a reasonable normaliza- 
tion. 

In this respect, looking for a more quantitative infor- 
mation and, as before, from the fit with g^ aa . = 2 GeV 2 
and based on the experimental procedure j37j . we have 
generated 8 differential cross section points with errors of 
1%, in the region 0.95 < q 2 < 1.3 GeV 2 and fitted the 
points with the exponential form, as shown in Fig. 12. In 
this case we have obtained 

A = 0.123 ± 0.004 mbGeV" 2 , 
B = 3.554 ±0.030 GeV" 2 , 

with x 2 /DOF = 2.52 for 6 degrees of freedom. We note 
that a close looking at the generated points in Fig. 12 
shows that the last four points have a slope slightly greater 
than the first four points and this effect seems also to be 
present in the measured elastic rates [3*6ll3"7] . Although the 
experimental data are still being analyzed by the DZero 
Collaboration, we understand that the above information 
and the numerical results displayed in Table 5, can con- 
tribute with the discussion on a suitable normalization for 
these elastic rates. We shall return to this point in what 
follows. 



4.3.2 Phenomenological aspects 

We now turn the discussion to the phenomenological con- 
text, with main focus on the result we have obtained 
for the highest energy with differential cross section data 
available, namely scattering at y/s = 1.80 TeV, Fig. 13. 
The point is to compare this result with predictions from 
the models referred to in the beginning of this Section. 

From Fig. 13, our result indicates a change of curvature 
in the region of the last three experimental points (q 2 « 
0.55 — 0.65 GeV 2 ) with a shoulder shape and not a dip 
(minimum) with defined position. This effect is due to the 
contribution from the real part of the amplitude as shown 
in Fig. 14, where we display separately the contributions 
to the differential cross section from only the real and 
only the imaginary parts of the amplitudes in the cases of 
Imax — 2 GeV 2 and q 2 nax — 14 GeV 2 . From this Figure 
we see that, as expected, the imaginary part presents a 
zero (change of sign) and inside this region, the value of 
the minimum in the differential cross section is due to the 
contribution of the real part (a shoulder in this case) . The 
real part of the amplitude also presents a zero at q 2 ss 
0.30 GeV 2 in the case of q 2 max = 2 GeV 2 and q\ « 0.38 
GeV 2 for g 2 ^ aa . = 14 GeV 2 . These results for the real part 
are in agreement with a theorem demonstrated by Martin, 
which states that the real part changes sign at q 2 > 0.1 
GeV 2 |38. 

On the other hand, the contributions from the imagi- 
nary parts are very similar in both cases, indicating a zero 
at ql w 0.70 GeV 2 for q 2 max = 2 GeV 2 and at q 2 « 0.73 
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Fig. 14. Contributions to the pp differential cross section at 
^J~s = 1.80 TeV from the real (solid) and imaginary (dashed) 
parts of the amplitude for q^ax = 2 GeV 2 and q^ax 14 GeV 2 . 



GeV 2 for q?^ ax = 14 GeV 2 . Therefore, from this Figure, 
we can infer with some security that the position of the 
first minimum in the differential cross section at this en- 
ergy occurs at ql = 0.70 GeV 2 (q 2 lax = 2 GeV 2 ). In the 
phenomenological context this value is in agreement with 
the predictions of the DGP, BSW and DL models, but 
not with that from the BGHP model, since the minimum 
in this model is predicted to be at g 2 «0.60 GeV 2 (coin- 
cident with the highest E710 point). We think that this 
is an important point, that should be carefully analyzed, 
when comparing elastic rates with model predictions at 
y/s = 1.96 TeV. 

Another aspect to note in Fig. 14 is that, in both cases, 
the contribution of the imaginary part dominates in the 
region of small momentum transfer, up to the beginning 
of the shoulder. On the other hand, in this region and for 
higher values of the momentum transfer, it is the contri- 
bution of the real part that dominates. However, in order 
to investigate this effect in more detail, we must consider 
the region of medium and large momentum transfer, that 
is the results of the fits with q^ax = 14 GeV 2 . We stress 
that, even under restrictive formal justification, our results 
taking into account all the differential cross section data 
are quite good, as shown in Fig. 9 and therefore it may 
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Fig. 15. Contributions to the differential cross section from 
the real (solid) and imaginary (dashed) parts of the amplitude 
for q 2 max = 14 GeV 2 . 



be instructive to discuss the implications of this variant of 
the fit. 

Concerning the contributions to the differential cross 
sections from the real and the imaginary parts of the am- 
plitude, we consider three typical examples: the results 
for pp scattering at 52.8 GeV and pp at 53 GeV, shown 
in Fig. 15, and those for pp at 546 GeV, displayed in Fig. 
16, together with the corresponding experimental data. 
The point is that, according to our predictions, in the en- 
ergy region of the CERN ISR (y/s w 23 - 63 GeV), the 
imaginary part dominates at medium and large values of 
the momentum transfer (Fig. 15). On the other hand, at 
higher energies, such as the regions of the CERN Collider 
(Fig. 16) and Tevatron (Fig. 14), it is the contribution 
from the real part that dominates. 

At last, it may also be instructive to see what kind 
of results can be predicted in the region of large momen- 
tum transfer at the RHIC, Tevatron and LHC energies. 
We display that in Fig. 17 up to q 2 = 8 GeV 2 , the in- 
terval generally considered in the publications. The main 
point here is the prediction of a smooth decrease of the 
differential cross section above the first minimum, with- 
out secondary structures in this region. Among the quoted 
phenomenological approaches, this behavior is predicted 
only in the DL model. However, we note from Fig. 9 that 
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a small change in the curvature is predicted at q 2 w 12 
GeV 2 . 



5 Conclusions and Final Remarks 

We have introduced an analytical paramctrization for the 
clastic hadron-hadron scattering amplitude and a fit pro- 
cedure characterized by at least five important novel as- 
pects: (1) the parametrization is almost model-independent, 
with enclosed dependences on the energy and momentum 
extracted from the empirical behavior of the experimen- 
tal data and in agreement with some high-energy theorems 
and bounds from AQFT; (2) the real and imaginary parts 
of the amplitude are entire functions of the logarithm of 
the energy s and are connected through derivative disper- 
sion relations; (3) the pp and pp scattering are also con- 
nected to the extent that analyticity and unitarity lead to 
dispersion relations; (4) the approach is predictive in both 
energy and momentum variables; (5) fits to pp and pp ex- 
perimental data, above 20 GeV, on the forward quantities 
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and then including differential cross sections, allow good 
global descriptions of all the data, even in different regions 
of the momentum transfer. 

We have presented a critical remark on a drawback 
that still remains in the present formulation, which is re- 
lated to the results for pP p (s), in the particular case of 
forward fit. One way to address this question may be to 
consider the derivative dispersion relations up to second 
or third order in the tangent operator. Results on this 
direction will be reported elsewhere. 

Another aspect that deserves some comment is the 
number of free parameters involved in the analysis. When 
including the differential cross section data, even in the 
region q 2 < 1 GeV 2 , the fit demands 3 exponentials in the 
imaginary part of the amplitude and therefore 30 free fit 
parameters. We understand that this can not be seen as an 
disadvantage of the formalism in terms of a large number 
of parameters. In fact, we are not treating with a theo- 
retical model but, on the contrary, a model- independent 
approach aimed to describe and predict the physical quan- 
tities of interest on empirical and formal grounds. There- 
fore, the number of parameters does not matter and, in 
this context, it can be as large as it is needed. 

In this analysis we made use of the standard sets of 
experimental data on pp and pp scattering above 20 GeV 
(referred to in Sect. 3.1), without any kind of data se- 
lection. As commented, this strategy explains the large 
values of the x 2 /DOF in the fits. However it is important 
to mention that recent analyses point out the necessity of 
some screening criterion in order to avoid spurious data, 
normalization problems and other effects in both forward 
and nonforward data .':!!) All that could improve the 
quality of the fits and will be subject of future investiga- 
tion. 

We also note that we did not use any model informa- 
tion in the construction of parametrizations (5-8): they 
were inferred only with basis on the empirical behavior 
of the experimental data above 20 GeV. However, from a 
phenomenological point of view, it is expected that some 
contributions from lower energies may still be present at 
the above threshold (for example, secondary mesonic ex- 
changes, in a Regge context @U). Therefore it may be in- 
teresting to test additional terms in our original parametriza- 
tion, that could simulate these effects, from an empirical 
point of view, and investigate the consequences in the de- 
scription of the experimental data. 

We now summarize some results we understand are 
topical in this analysis. The behavior of the forward quan- 
tities, <Jtot(s) and p(s), from pp and pp scattering are char- 
acterized by crossing effects, which are typical of Odderon 
contributions. A relevant result is the prediction that p pp 
becomes higher than fP v above yfs w 80 GeV, a result 
that might be verified in short term at the RHIC by the 
pp2pp Collaboration. Our results for the differential cross 
section at the Tevatron energies are in agreement with the 
predictions from the majority of the models, except that 
by Block et al. [35], in what concerns the position of the 
first minimum. We have also discussed the applicability of 
our numerical results in the normalization of elastic rates. 



We add that, if we consider the fit including all the differ- 
ential cross section data, the DL model is favored, since no 
structures is predicted in the region of large momentum 
transfer. 

In closing we should stress that, despite the encourag- 
ing results we have reached, this phenomenological analy- 
sis constitutes a first attempt in the search of a formally 
rigorous and predictive model-independent approach. Much 
more research must still be done along several lines, as for 
example, a complete check on all the high-energy theo- 
rems and bounds, to establish the exact interval in the 
momentum transfer variable in which dispersion relations 
hold (or another framework for evaluation of the real part 
of the amplitude), studies on the effect of higher orders in 
the derivative dispersion relations and a systematic inves- 
tigation on the influence of data selection. We hope that 
the results here presented can contribute with further de- 
velopments along these aspects. 
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